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The phase transition from a superfluid to a Mott insulator has been observed in a 23 Na Bose- 
Einstein condensate. A dye laser detuned ~ 5 nm red of the Na 3 2 S^ 3 2 Pi/2 transition was used 
to form the three dimensional optical lattice. The heating effects of the small detuning as well 
as the three-body decay processes constrained the timescale of the experiment. Certain lattice 
detunings were found to induce a large loss of atoms. These loss features were shown to be due to 
photoassociation of atoms to vibrational levels in the Na2 (1) 3 X1^ state. 

PACS numbers: PACS 03.75.Lm, 73.43.Nq, 33.20.Kf, 34.20.Cf 



Introduction 



Optical lattices have become a powerful tool to en- 
hance the effects of interaction in ultracold atomic sys- 
tems to create strong correlations and probe many-body 
physics beyond the mean- field theory 0, 0, IE EE IE IE E_f • 
Simply through varying the depth of the lattice poten- 
tial, one changes the tunneling rate as well as the on-site 
interaction energy by changing the confinement of the 
atoms. The strength of the atomic interaction can be di- 
rectly tuned with a magnetic Feshbach resonance In 
comparison to 87 Rb, which has been used in almost all 
experiments on optical lattices, 23 Na has stronger and 
wider Feshbach resonances that are experimentally ac- 
cessible |E [H3- ne such resonance has been used to 
produce quantum degenerate Na2 molecules [l]| • There- 
fore, a sodium condensate loaded into an optical lattice 
would be a rich and flexible system for studying strong 
correlations. 

So far, most optical lattice experiments have been per- 
formed with relatively heavy atomic species (e.g. ru- 
bidium and potassium) for which the recoil frequencies 
are lower and lasers are readily available to achieve trap 
depths of several tens of recoil frequencies at a few tens 
of milliwatts. For 23 Na, high power single-mode lasers 
are necessary for similar experiments. In this work, we 
chose to use a dye laser red-detuned by ~ 5 nanometers 
from the D lines of sodium (589 nm). The spontaneous 
scattering rate limited the time window of the experi- 
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ment to less than 50 ms, but was still sufficient to satisfy 
the adiabaticity condition to explore the quantum phase 
transition from a superfluid to a Mott insulator. We also 
observed strong atom losses at various lattice laser de- 
tunings, which were interpreted as photoassociation tran- 
sitions. The particular molecular states responsible for 
these transitions were identified through theoretical cal- 
culations and previous experimental data. 



Experiment Setup 

A 23 Na Bose-Einstein condensate containing up to 10 6 
atoms in the \F = 1, = —1) state was first produced 
in a magnetic trap and subsequently loaded into a crossed 
optical dipole trap. The optical trap was derived from 
a single-mode 1064 nm infrared laser, with the hori- 
zontal and vertical beams detuned by 60 MHz through 
acousto-optic modulators. The number of condensed 
atoms was varied through three-body decay in a tight 
trap (v x ,y,z = 2tt x 200,328,260 Hz) , after which the 
trap was decompressed (oJ x ,y,z = 27r x 110, 155, 110 Hz) 
to allow further evaporation and re-thermalization. A 
vertical magnetic field gradient was applied to compen- 
sate for gravity and avoid sagging in the weaker trap. 

A dye laser operated at 594.710 nm was used to set up 
a three dimensional optical lattice. The three beams were 
focused to l/e 2 -waist of ~ 82 /im at the condensate, and 
retro-reflected to form standing waves. The two horizon- 
tal beams were orthogonal to each other, while the third 
beam was slanted at ~ 20° with respect to the vertical 
axis due to limited optical access. The three beams were 
frequency-shifted by ±30 MHz and 80 MHz to eliminate 
cross interference between different beams. 



2 



The gaussian profile of the lattice beams added an ad- 
ditional harmonic trapping potential, while the localiza- 
tion of atoms at the lattice sites increased the repulsive 
mean field interaction. At the maximum lattice depth, 
the trap frequencies due to the combined potentials of the 
optical dipole trap and the lattice beams were ~ 510 Hz 
for all three dimensions. The trap parameters were cho- 
sen such that during the ramping of the optical lattice 
potential, the overall size of the cloud (parametrized by 
Thomas-Fermi radii) remained approximately constant 
in order to minimize intra-band excitations (the mean 
Thomas- Fermi radius is ~ 14 jam for 10 6 atoms). The 
peak per-lattice-site occupancy numbers achieved in our 
experiment were between 3 to 5. 



Quantum Phase Transition 

Atoms held in a shallow optical lattice can tunnel freely 
from site to site and form a superfluid phase. As the lat- 
tice is made deeper, the atomic interaction is increased 
while the tunneling rate between lattice sites is exponen- 
tially suppressed. The system then undergoes a phase 
transition to an insulating phase - the Mott-insulator - 
in which each lattice site contains a definite fixed num- 
ber of atoms. According to the mean- field theory for the 
homogenous systems of atoms in the lowest band of an 
optical lattice, the critical point for the phase transition 
from a superfluid to a Mott-insulator state with n atoms 
per lattice site is determined by [13 . H3 . : 



where: 



U = z(2n + 1 + 2y/n(n + l))J 



U = / d s x\w(x) 4 

m 



(i) 



(2) 



is the on-site interaction energy; 
J = Jd 3 xw*(x- Xiatt/2)(-^V 2 + V latt (x))w(x) (3) 

is the tunneling rate between adjacent lattice sites; z is 
the number of nearest neighbors in the lattice (6 for a 
cubic lattice); m is the atomic mass; a s is the s-wave 
scattering length (2.75 nm for 23 Na)); w(x) is the Wan- 
nier function; Xi a tt is the lattice wavelength; Vi a tt{x) is 
the lattice potential. 

Figure[I]shows U and J n = z(2n^l^2^/n(rl^T))J for 
a cubic lattice as a function of the lattice depth, obtained 
through a band-structure calculation. All energies are ex- 
pressed in units of the recoil energy E reco u = h 2 kf att /2m, 
where ki a tt — ^ I Matt is the lattice wavenumber. With 
this scaling J is independent of Matt- The peak occu- 
pancy number in our experiment was < 5. From Fig. 
we find that the the critical points are at a lattice depth 
of 14.2, 16.2, 17.6, 18.7, and 19.5 (all in units of E recoU ) 



for n = 1, 2, 3, 4, and 5 respectively. The inset of Fig. [I] 
shows that the ratio of U/J increases rapidly with in- 
creasing lattice depth. 

When a weak harmonic trap is present in addition to 
the lattice potential, as is the case for the experiment, 
the atomic density is not uniform. Nevertheless, Eqs. (0 
OJ) can be used to estimate the lattice depth needed to 
observe the Mott-insulator phase transition at any point 
in the harmonic trap. Given the local density of the 
initial condensate, a local value of n can be estimated 
and thus the local critical lattice depth can be read off 
from Fig. Since the critical depth increases with n, 
one expects that as the lattice depth is increased, shells 
of different occupancies will undergo the transition to 
the Mott-insulator phase starting from the edge of the 
density profile and moving in towards the center. 
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FIG. 1: Determination of the phase transition points: The 
phase transiti on at occ upancy number n occurs when J n = 
z(2n + 1 + 2^/n(n + 1)) J as given by equation (J equals U. 
The horizontal location of the crossing point where J n = U 
gives the critical lattice depth. The inset shows the ratio of 
U over J as a function of the lattice depth. 

In our experiment, the optical lattice was linearly 
ramped up to a maximum potential of 20 E recoi i in a 
variable time r rarnp (E recoi i = h2AA kHz for our sys- 
tem). The lattice depth was calibrated by probing the 
energy difference between the first and the third band 
at zero quasi-momentum with small amplitude modula- 
tion of the lattice beams (see, e.g., 6]). After reaching 
the peak value, the lattice was ramped back down again 
in r rarnp . The ramp sequence was stopped at different 
times when both the trap and the lattice were suddenly 
switched off (in < 1 /is). Absorption images were then 
taken after some time-of- flight as shown in Figure O The 
disappearance of the interference pattern as the lattice 
depth was increased indicated the loss of phase coherence 
and a transition from the superfluid state to the Mott 
insulator state P|. The subsequent revival of the inter- 
ference patterns as the lattice depth was reduced ensured 
that the system remained essentially in the ground state 
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during the ramping process. Different r rarnp s were used 
to check the adiabaticity condition. The peak sponta- 
neous light scattering rate was about 21 s _1 at the max- 
imum intensity. Therefore for r rarnp < 10 ms, less than 
20% of the atoms spontaneously scattered a photon. 

After the lattice was fully ramped down, most of the 
atoms (> 80%) remained in the condensed fraction while 
the rest were heated and distributed across the first Bril- 
louin zone. Based on the number of atoms that remained 
in the condensate after the lattice was fully ramped down, 
we conclude that r rarnp > 1 ms satisfies the intra-band 
adiabaticity condition. In the following discussion, all 
measurements were performed for r rarnp = 1,5,10 ms, 
but only the data for r rarnp = 5 ms are shown as repre- 
sentative of similar results unless otherwise noted. 

To characterize the lifetime of the Mott-insulator 
phase, we held the lattice depth at the maximum level 
for various amounts of time before ramping the lattice 
down to 8 E reco u (below the Mott-insulator transition 
point) and taking the time-of-flight image. If the system 
remains in the ground state, the contrast of the inter- 
ference pattern should be recovered, whereas additional 
heating populates the Brillouin zone and reduces the in- 
terference contrast. A cross-section of the density profile 
was taken along the horizontal direction showing the in- 
terference peaks on top of a broad background (see Figure 
0J). The 5 interference peaks and the broad background 
were fit by 6 gaussians. The ratio between the total inte- 
grated area of the peaks and the background was used as 
the contrast to quantify the heating of the system. The 
contrast gives a more sensitive measure of the heating 
compared to simply counting the recovered condensate 
atoms. As the atoms in the interference peaks quickly 
move apart, they are not as broadened by the mean- field 
expansion as a single condensate. 

We performed the same measurement for two different 
peak occupancy numbers n ~ 3 and 5. FigureOlshows the 
decay of the contrast and the lifetime r was determined 
using an exponential fit. The fitting error on the lifetime 
was less than 17%. The lifetime was about 50% longer 
for n = 3, implying that inelastic collision processes sig- 
nificantly contributed to the heating of the system. The 
three-body decay rate at the maximum lattice potential 
for the peak on-site atomic density (~ 10 16 cm -3 for 
n = 5) is about 100 s _1 0, consistent with the observed 
lifetimes of ~ 10 ms. The peak density of a condensate 
in a harmonic trap and therefore the peak occupancy 
number scales with 2/5 power of the total number of 
atoms, and our method for varying the number of atoms 
(through three body decay) was unable to consistently 
produce low enough atom numbers for peak occupancy 
< 2. The signal-to- noise ratio of our current imaging sys- 
tem also became marginal for such low atom numbers. 



Photoassociation Resonances 

In this experiment, in addition to losses due to three 
body recombination, we observed large losses of atoms for 
certain specific tunings of the lattice laser in the range 
592 nm to 595 nm. A sample of such a loss feature is 
shown in Fig. For this measurement, the same ramp 
sequence was used as before with r rarnp = 1 ms. The 
peak intensity is about 280 W/cm 2 in each lattice beam. 
Due to the intentional frequency shifts between the three 
lattice beams the effective bandwidth of the lattice light 
as seen by the atoms is ~ 100 MHz. For the narrow fre- 
quency scan range of Fig. 0] the relative frequency scale 
was determined to better than 25 MHz using a Fabry- 
Perot cavity with a 2 GHz free spectral range. 

Single-photon photoassociations proved to have caused 
these losses. The lattice laser is tuned by ~ 160 cm -1 to 
the red of the atomic 3 2 S^ 3 2 P3/2 transition and thus 
is in a spectral region where it might drive photoasso- 
ciation transitions [2, El to rovibrational levels in 
molecular states dissociating to either the 3 2 S+3 2 Pi/ 2 or 
3 2 S+3 2 P3/2 limits. Such a photoassociation transition, 
followed by the spontaneous radiative decay of the ex- 
cited molecule into either a bound ground electronic-state 
molecule or into "hot" atoms, results in significant losses 
of atoms from the lattice. It is therefore important to 
identify the locations and strengths of these resonances 
and choose the appropriate lattice wavelength to avoid 
such losses. 

There is an extensive body of knowledge on the pho- 
toassociation of ultracold alkali-metal atoms and the be- 
havior of the molecular potentials dissociating to the 
3 2 S+3 2 P 1/2 or 3 2 S+3 2 P 3/2 limits 000. Figure 
shows the relevant excited molecular potentials as a func- 
tion of internuclear separation R. The ground electronic 
states of Na2 are the and a 3 E+ states and two col- 

liding ground state atoms will be some mixture of these 
symmetries. To the extent that the excited states are well 
described as S or II states, the g <-> u and AS — se- 
lection rules imply that photoassociation transitions are 
allowed only to the two E states (A 1 Ej and (1) 3 E+). 

Previous experiments have identified the locations of 
the strong transitions to the A X E+ stat e I l9l and the 
weak transitions to the (l^Hgfig) state |2fl|. These are 
shown in Fig.[BJ We looked for but failed to find any sig- 
nificant losses attributable to the weak (l) 1 n^(l^) state 
resonances. We were able to confirm one of the A- state 
resonances, indicated by the dot in Fig. Since the 
primary goal of the present work was to avoid photoas- 
sociation losses we did not investigate known A-state lo- 
cations further. Our search focused on the strong res- 
onances due to the (1) 3 E+ state which had not been 
previously observed in this spectral region. In 0.5 GHz 
steps, we scanned through a 30 GHz range around the 
theoretically predicted locations based on the model and 
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FIG. 2: (a) Observation of the superfluid to Mott insulator transition: The lattice depths for the sequence of images are from 
left to right 0, 4, 8, 12, 16, 20, 16, 12, 8, 4, E recoi i. The time-of-flight is 7 ms and the field of view is 1000 /mix 1200 /mi. The 
peak occupancy number n is 5 for these images, and the phase transition should occur between 14 and 20 E reC oii according to 
the mean-field calculation, (b) Time dependence of the lattice depth. 



auxiliary experimental data discussed below. In all but 
one such scans, we were able to observe between 1 to 3 
dips in the remaining atom numbers within a range of 
< 15 GHz, including the loss feature shown in Fig. 21 
As shown in Fig. the agreement of the observed loca- 
tions with the predictions and auxiliary measurements 
confirms that these loss features are due to photoasso- 
ciation to the (1) 3 E+ state. The locations of the rovi- 
brational levels of the b 3 II M state are not known in the 
current tuning range. Their spacing, however, should be 
equal to that of the (l) 1 !^ levels and our observations 
are not consistent with such spacings. 

The potential curves used in the calculation of the 
(1) 3 H+ vibrational levels were generated from an ex- 
tended version of the model developed by Movre and 
Pichler for calculating the combined effects of the 1/R 3 
resonant dipole interaction and the atomic spin-orbit in- 
teraction. Such models have been extensi vely used to 
interpret photoassociation experiments [HJ 0, [H| • To 
the long range potentials generated by the Movre- Pichler 
type calculation we append the results of ah initio cal- 
culations on the short range molecular potentials (in the 
chemical bonding region). These short range potentials 
are not sufficiently accurate to allow predictions of ab- 
solute vibrational positions. It is necessary to make 
slight adjustments to these potentials to match experi- 



mentally measured vibrational positions, which were ob- 
tained through a separate photoassociation spectroscopy 
experiment in a dark-spot magnetooptical trap (MOT) 
containing Na 3 2 S(/ = 1) atoms at about 300 /iK, us- 
ing a two color ionization scheme and an ion detector 
22, 23]. The photoassociation spectra taken in the MOT 
have higher resolution than the loss features in the lattice 
experiment. 

In the spectral region of interest for the lattice exper- 
iment, three vibrational levels of the (1) 3 £+ were iden- 
tified by measurements in the MOT. The spectra show 
a 1 9 component with a complicated hyperfine/rotational 
pattern, and, slightly higher (~ 0.3 cm -1 ) in energy, a 
0~ component with a simpler nearly rotational pattern. 
This ordering of the l g and 0^ components is in agree- 
ment with the Movre-Pichler model. For these levels, 
the J = feature of the 0~ component were found to 
be at 16845.155 cm" 1 , 16852.585 cm" 1 , and 168908.201 
cm -1 , calibrated to iodine lines with an estimated un- 
certainty of about 0.004 cm -1 and are shown in Fig. 
This new data on the (1) 3 S+ was sufficient to calibrate 
our extended Movre-Pichler model without any further 
adjustable parameters. The predicted positions shown 
in Fig. agree well with those of the loss features ob- 
served in the lattice experiment, thus identifying the loss 
features as photoassociation transitions to levels in the 
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FIG. 3: Lifetime of the Mott insulator state: The top row 
of absorption images were taken for hold times 0, 5, 10, 15, 
20 ms. The second row is the horizontal cross section of the 
atomic density profile taken at the center of the absorption 
images. The decay of the Mott insulator state is observed as 
the loss of contrast of the interference peaks. The curves are 
exponential fit. The error bars in the graph are statistical 
standard deviation of multiple shots. 

(1) 3 S+ state. 

Based on this insight, we chose a wavelength of 
594.710 nm for our lattice experiment (corresponding to 
-158.5 cm -1 in Fig. EJ. This tuning lies > 45 GHz from 
the closest molecular resonance. Given the observed on- 
resonance photoassociation rate of 1 ms -1 in the lattice, 
photoassociative decay can be ignored at such detunings 
as the rate scales as the square of the ratio of the natural 
linewidth to the detuning. 



Conclusions 

In this paper, we explored the possibility of using a 
dye laser detuned « 5 nm from the Na D lines to study 
many-body physics of a sodium BEC in an optical lat- 
tice, which could allow for an independent control of 
interaction using magnetic Feshbach resonances. The 
superfluid to Mott-insulator transition was observed in 
a Na Bose-Einstein condensate for the first time. The 
main technical difficulties are due to the heating from the 
spontaneous light scattering and three-body decay pro- 
cesses. In addition, several photoassociation resonances 
were observed and identified by means of auxiliary spec- 
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FIG. 4: Photoassociation resonances: Atom loss from the op- 
tical lattice as a function of the lattice laser detuning. The 
top row of images shows the remaining atom cloud after the 
lattice was turned off and the cloud expanded for 30 ms. The 
lower graph shows the remaining number of atoms for these 
and additional frequency points. The zero of the x axis cor- 
responds to a laser wavelength of 594.490 nm as determined 
with a wavemeter with accuracy ~ 1 GHz. The error bars are 
one standard deviation uncertainty. The line is only a guide 
and does not have a theoretical underpinning. 



troscopy measurements combined with theoretical mod- 
eling. These resonances were avoided by choosing an ap- 
propriate lattice wavelength. In future experiments, we 
plan to use a high-power single-mode infrared laser at 
1064 nm to eliminate atomic light scattering (Figure Q] 
also shows the transition points for a 1064 nm lattice). 
Moreover, heating from three body recombination can be 
avoided by using occupancy numbers less than three. 
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FIG. 6: Comparison of loss features (plusses and dot) observed in a Bose condensate trapped in an optical lattice, to known 
molecular level positions: The x axis is the binding energy of the molecular level relative to the 3 2 S(/ = 2) + 3 2 P 3 / 2 (/ = 3) 
atomic asymptote. If a loss feature, observed at a laser photon energy of ^photon, is due to photoassociation of two 3 2 S(/ = 1) 
Na atoms to a bound molecular level, then the molecular level is bound by E — ^photon where ^o/(^c)=16973.4636 cm -1 . 
The laser tunings of observed loss features are plotted in this way (the pluses and the dot). The lines labeled (l) 1 n p (l p ) are 
binding energies of the J — 1 rotational levels of vibrational levels in this potential measured in a magnetooptical trap (MOT) 
f |2(1| and unpublished work). The lines labeled A 1 S^(0j) are the energies reported in the RKR curve of (nominally 
J — positions) shifted to agree with the photoassociation measurement of 21] corrected for rotation. The lines labeled 
(l) 3 Xj~ are from an extended Movre-Pichler model calculation calibrated to experimental data (the crosses) from a two-color 
photoassociation experiment in a MOT. The lines displaced upwards (downwards) are the l g (0~) component of the (1) 3 E+ 
state. The crosses indicate the three calibration points which lie in this spectral range and are measurements of the J — 
levels of the 0~ component. 



